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In equilibrium, at finite temperature below and above the deconfining phase transition, we have 
generated lattice SU{2) gauge fields and have exposed them to smearing in order to investigate the 
emerging clusters of topological charge. Analysing in addition the monopole clusters according to 
the maximally Abelian gauge, we have been able to characterize part of the topological clusters to 
correspond either to non-static calorons or static dyons in the context of Kraan-van Baal caloron 
solutions with non-trivial holonomy. We show that the relative abundance of these calorons and 
dyons is changing with temperature and offer an interpretation as dissociation of calorons into dyons 
with increasing temperature. The profile of the Polyakov loop inside the topological clusters and the 
(model-dependent) accumulated topological cluster charges support this interpretation. Above the 
deconfining phase transition light dyons (according to Kraan-van Baal caloron solutions with almost 
trivial holonomy) become the most abundant topological objects. They are presumably responsible 
for the magnetic confinement in the deconfined phase. 

PACS numbers: 11.15.Ha, ll.lO.Wx 
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I. INTRODUCTION 



In our previous work P| we have started an investi- 
gation of topological objects that appear in the course 
of smearing ^ of equilibrium lattice fields. These were 
generated at finite temperature in the confining phase of 
SU(2) lattice gauge theory. The main idea was to analyze 
the monopole content of these objects, still far from being 
classical solutions, according to the maximally Abelian 
gauge (MAG) in order to classify them as related to 
non-static calorons or static dyons. Such a classification 
would be natural in the context of Kraan-van Baal (KvB) 
caloron solutions with non-trivial holonomy &| . In 

the present paper, using some additional techniques, we 
take up this problem again, now for a couple of temper- 
atures below and above Tdoc- 



Before pointing out the new tools for this analysis, we 
should briefly recall the general view and the terminol- 
ogy used, that originate from the classical KvB solu- 
tions 0, 01 • We will do this for the simplest case, the 
SU{2) gauge theory. In this case, at finite temperature, 
the classical carriers of one unit of topological charge 
are selfdual or antiselfdual objects (periodic caloron so- 
lutions) which consist of two monopole "constituents" . 
In the limit of small inter-monopole distances they form 
a single, non-static caloron, whereas at large distances 
they dissociate into a pair of separate static lumps of ac- 
tion and topological density (static Bogomolnyi-Prasad- 
Sommerfield (BPS) monopoles or "dyons"). In the 
latter case the action Sn'^/g'^ (or one unit of topological 
charge) is shared among the constituents according to 
the holonomy of the gauge field, in the ratio ZJ/uj, where 



uj = 1/2 — cj. The asymptotic value of the Polyakov loop 

Loo = Jim L{x) 

\x\ — *oo 

= lim - tr P exp | i j A4{x,t) dt ) 

\x\ — >oo 2 







cos{2ttlu), 0<cj<1/2, 



(1) 



determines the holonomy parameter to and its comple- 
ment ZU. 

Each physical phase can be thought as a medium creat- 
ing a certain Loo as boundary condition for the calorons. 
For example, maximally non-trivial asymptotic holon- 
omy. Loo = or UJ — oj = 1/4:, forces the constituents to 
carry equal topological charge ±1/2. With the holonomy 
becoming more and more trivial with increasing temper- 
ature, Loo ^ il , the constituents become imbalanced. 
On the other hand, the constituents are distinguished 
by the local behavior of the Polyakov loop L{x) w ±1 
close to their positions. This is the characteristic signa- 
ture, irrespectively whether the action forms two sepa- 
rated lumps or a single lump, the latter one with unit 
topological charge. 

Viewed by Abelian monopoles appearing in the MAG, 
static dyons are represented by monopoles temporally 
wrapping around the lattice, whereas single non-static 
calorons are localized in time and have monopole loops 
of finite extent running close to the center. This feature 
of calorons in the MAG has been pointed out first by 
Brower et al. . This was our guiding hypothesis already 
m Rcf. "ll. In Fig. n we show a classical caloron solution 
in the two limiting cases together with the corresponding 
MAG monopole content. In Ref. ^ we have seen such 
static monopoles also in the Polyakov gauge, regularly 
accompanying static, i.e. dissociated caloron configura- 
tions while being rarely found in non-static calorons. 

According to the scenario to be outlined in Sec. HTl we 



FIG. 1: (color online) Left: two static dyons accompanied by two static monopoles; right: a single caloron accompanied by a 
monopole loop. The classical SU{2) configurations have been generated on a 16^ x 4 lattice and then Abelian-projected. In the 
{z,y,t) lattice at x = we show the sites with action density s > Smax/5 (left, for dyons) and s > Smax/40 (right, for calorons) 
together with the emerging monopole trajectories. 



expect to observe below Tdcc both dissociated and non- 
dissociated clusters of topological charge ±i and ±1, re- 
spectively, in some mixture. Above Tdoc only clusters 
with charge ±1 should remain, with a density decreasing 
with increasing temperature. That means that we are 
expecting to see calorons with trivial holonomy actually 
dominating the topological charge above Tdcc- 

The topological structure in the deconfined phase is 
partly also a question concerning the capability of the 
smearing method applied in this paper. Whereas the 
topological structure is known to be preserved under 
smearing in the confined phase, experience shows that 
under comparable smearing conditions the topology is 
rapidly wiped out in the deconfined phase. Therefore, in 
order to escape these bad prospects, we have to modify 
the smearing conditions in the deconfined phase. 

We point out in Sec.|n]that there are fermionic meth- 
ods to determine the topological density but, in order 
to detect a semiclassical structure, certain UV smear- 
ing techniques have to be applied as well. Our analo- 
gous method applied for the same purpose is a combina- 
tion of gluonic measurements of the most naive topologi- 
cal charge density combined with the monopole localiza- 
tion in the maximally Abelian gauge. Both are applied 
to APE smeared configurations. In a very recent pa- 
per it has been demonstrated, by comparison of smeared 
and cooled configurations with the original equilibrium 
configurations, that the spectrum of a chirally improved 
Dirac operator remains practically unaltered meaning 
that the infrared structure "seen by the low-lying chiral 
fermion modes" in equilibrium, is conserved under con- 
trolled cooling or moderate smearing. 

Since for our purpose the correlation with Abelian 
monopoles is of crucial importance, the technical im- 
provements (compared to Ref. 1]) in the present pa- 
per are concerning the methods detecting the monopole 
structure. This applies both to gauge fixing and to the 



analysis of separable monopole loops. We have now fixed 
the MAG by the more advanced technique of simulated 
annealing [lOj. This method is known to lead, in compar- 
ison with the overrelaxation method, to higher maxima 
of the gauge functional and to a minimal monopole con- 
tent. Analogously to the preceding paper, the detection 
of monopole trajectories is used here to identify clusters 
of the topological charge density resembling the two ex- 
treme appearances of KvB calorons. Unlike the previous 
paper, we are now actually exploring the full monopole 
content of individual topological objects, not only look- 
ing for intersections of time-like monopole currents with 
the topological charge clusters 0. We are taking the 
precise type of monopole clusters into account. 



Last but not least, we use now a larger lattice (with size 
24^ X 6) and a certain set of different /3- values. Three val- 
ues, P — 2.2, 2.3, 2.4, belong to the confined phase while 
the remaining two values oi P — 2.5, 2.6, represent the 
deconfined phase. This will allow us to get information 
about the temperature dependence of the caloron/dyon 
composition that was left open in the previous investiag- 
tions. 



In Sec. ^we will outline the general physical picture 
behind the present investigation. In Sec. IIIII we will 
briefiy explain the set of our runs and in particular the 
conditions for the smearing procedure. In Sec. IIVI we 
shall describe the local infiuence of monopoles altering 
the distribution of the Polyakov loop. This part is in- 
teresting by its own and is not related to the following 
cluster analysis. In Sec. 0we will describe the results 
of the analysis of topological clusters. The classification 
with the help of monopoles will be of use here. We draw 
our conclusions in Sec. IVII 
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II. THE PHYSICAL PICTURE 

The one-loop calculation of the caloron amplitude 
has shown that below some temperature calorons may 
become instable with respect to the separation into their 
constituents. It would be attractive to identify this tem- 
perature with the deconfinement temperature because, 
in a dilute gas calculation, it has been shown that at 
this temperature trivial holonomy Loo ~ 1 turns from 
a minimum of the free energy P (as a function of Loo) 
into a maximum. The global view of the function P{Loo) 
is, however, beyond the capability of the dilute caloron 
gas approximation. In a schematic dyon gas model, 
Diakonov 12] has demonstrated how the stability of 
Loo — (being the minimum of P) could eventually be 
achieved. 

All this taken together leads to a simplified picture 
that describes the essence of the transition to confine- 
ment as the dissociation of tightly bound calorons (with 
trivial asymptotic holonomy) into dissociated ones with 
nontrivial holonomy and discernible dyonic constituentis. 
This cannot be the full truth, however. As we know from 
Refs. 0, 13 J at low temperature calorons, even of non- 
trivial holonomy, are likely to appear as single lumps of 
action that nonetheless are different from the well-known 
instanton solutions because of the nontrivial holonomy 
boundary conditions they must satify in the confining 
medium. 

Therefore we expect to observe in a study at different 
temperatures the constituents becoming manifest inside 
dissociated calorons more abundantly approaching the 
deconfining transition from below. Whereas the topolog- 
ical susceptibility is temperature independent through- 
out the confined phase, the relative abundance of static 
objects with an appreciable fractional topological charge 
(close to ±i) should increase towards the phase tran- 
sition temperature compared to the abundance of non- 
static calorons. 

We have announced in the Introduction that we have 
used four-dimensional smearing for this investigation. 
There is an alternative way to locate calorons and their 
constituents, respectively. It is based on the localiza- 
tion behavior of the zero modes and the near-zero modes 
of the Dirac operator 0, A jumping behavior of 

the zero mode in charge Q = ±1 configurations can be 
observed under the influence of changing boundary con- 
ditions, deliberately imposed on the fermion field. This 
is very reminiscent to what is familiar from the case of 
an exact caloron background field. Interpolating with an 
appropriate phase factor between periodic boundary con- 
ditions and antiperiodic ones, the zero mode may jump 
from one constituent to another. 

For various classical solutions extracted on the lat- 
tice by the cooling method this has been established for 
SU{2) [13 as well as for SU[2,) For equihbrium 

lattice gauge fields selected to possess unit topological 
charge, in the SU{'i) and SU{2) cases a similar "jump- 
ing" of lowest-lying modes of a chirally improved lattice 



Dirac operator was reported by Gattringer and coau- 
thors Of course, Q = ±1 equilibrium con- 
figurations generically contain more clusters and spatial 
fluctuations of topological charge than the single zero- 
mode is able to detect. Such configurations differ from a 
pure (anti-)caloron even in the deconfincd phase. It has 
been verified '2^ , that the single zero mode is always at- 
tracted to positions that are maxima of the topological 
density inside extended topological clusters. This check 
has required a very well-tuned smearing and has used a 
topological charge operator defined in terms of the naive 
twisted plaquette. However, there were more clusters 
and cluster maxima going undetected by the single zero 
mode. In our previous 0] and in the present paper, the 
diagnostic role is shifted to the MAG monopoles. Note 
that these are not limited in number or to a particular 
chirality as the zero modes are. 

With a chirally perfect or improved Dirac operator at 
hand (that satisfies the Ginsparg- Wilson equation) , there 
is a definition [2l| of the topological charge density which 
is equivalent, in the continuum limit, to the gluonic def- 
inition. It may supersede a gluonic definition at finite 
lattice spacing. This density also exposes a clustering 
pattern |23 . l28[ |23| strongly dependent on ultraviolet fil- 
tering, usually implemented by a cutoff (in eigenvalues) 
of the contributing eigenmodes. The non-filtered den- 
sity, in contrast, contains fluctuations of all scales and 
cannot be related to a semiclassical model. An analysis 
of calorons and dyons using these tools along the ideas 
of the present paper is left for the future. 

In Ref. it has been attempted to construct single- 
caloron and multi-caloron solutions at lower temperature 
starting from high temperature ones. It was character- 
istic that all classical constituents were delocalized to a 
scale set exclusively by the periodic 4- volume, in contrast 
to the (jumping) zero modes seen also in zero tempera- 
ture equilibrium configuration background configurations 
which are highly localized due to quantum fluctuations 
of the gauge field. This limits our expectations to find, in 
the limit of vanishing temperature, a semiclassical struc- 
ture that would reveal the elusive "instanton quarks" in 
separation. 

Using the new procedures under wider physical con- 
ditions at finite temperatures, similar to Ref. Q we will 
observe both non-static calorons and static dyons. In the 
confined phase the relative population of isolated dyons 
will be seen growing with increasing temperature. This 
will be interpreted as progressive dissociation (dipole 
splitting) of calorons into dyons in the confined phase. 
This is the main result of our paper directly related to 
the KvB caloron picture and confirming the simplified 
picture described above. 

Unexpectedly from the point of view of the scenario 
outlined above, in the deconfined phase nondissociated 
calorons are rare objects, even though charge Q = 
±1 configurations are still observed after a well-tuned 
amount of smearing. According to the different asymp- 
totic holonomy (changing away from vanishing Polyakov 
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loop toward Loo 0) one would expect an extreme asym- 
metry among the constituent dyons if they appear at all. 
One type of dyons should be small in size and heavy in 
action with a peak of the Polyakov loop of opposite sign 
with respect to Lao, and another type of dyons large in 
size and light in action with a peak of the Polyakov loop 
of equal sign with the overall holonomy. The heavy dyons 
should be as rare as the nondissociated calorons are, but 
they are not found paired with the light dyons. The 
latter will be found about an order of magnitude more 
frequently and contribute less to the topological charge. 
Being static and always connected with a static Abelian 
monopole, they could be held responsible (together with 
other monopoles) for the magnetic confinement in the de- 
confined phase. This is the second main observation to 
be reported. 

III. DETAILS OF THE SIMULATIONS AND OF 
SMEARING 

We have generated 200 independent Monte Carlo con- 
figurations for each /? = 2.2, 2.3, 2.4, 2.5, 2.6 on a 24^ x 6 
lattice. We used the heatbath method according to the 
Wilson plaquette action. These samples characterize fi- 
nite temperature both in the confined and deconfined 
phases. At Nr — 6, the critical (3 corresponding to the 
deconfining phase transition is about (3dec = 2.42. For 
each (3 the corresponding ratio T/Tdoc is shown in Table 
I. The ensembles further have been subjected to smooth- 
ing by four-dimensional APE smearing 0- The fixed 
smearing parameter was a ~ 0.45, whereas various num- 
bers N of iterations were applied at different (3. From 
our previous observations on a 20'^ x 6 lattice at (3 — 2.3 
we concluded that dyons become visible above the noisy 
background after = 50 smearing steos. 
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FIG. 2: (color online) The effect of smearing on the "string 
tension" a. In confinement (/3 = 2.23, 2.3, 2.37) a was ob- 
tained from the correlator of Polyakov loops. In deconfine- 
ment (/? = 2.5, 2.6) a was obtained from the behavior of 
space-like Wilson loops. The lattice size was 20^ X 6. 



For the present investigation at /3 = 2.3 and the two 
other f3 = 2.2, 2.4 in the confined phase we adopt 
this number of smearing steps, justified by the obser- 
vation that the string tension ^ at the neighbouring val- 
ues f3 — 2.23, 2.3, 2.37 is reduced after 50 smearing 
steps to approximately the same 60% of the original value 
(see Fig. EJ. In the deconfined phase, where the spatial 
string tension was used for this comparison, the same per- 
centual reduction of this string tension has been observed 
at A^ = 25 and A^ = 20 smearing steps for (3 = 2.5, 2.6, 
respectively. 



IV. LOCAL POLYAKOV LOOP, MONOPOLES 
AND ASYMPTOTIC HOLONOMY 

Our first observation, as in Ref. (Ij, is the correla- 
tion between the values of the Polyakov loop and the 
presence of Abelian monopoles. For this purpose we 
have averaged the Polyakov loop over all eight corners 
of a three-dimensional cube (dual to a time-like link of 
the dual lattice) where an Abelian magnetic charge has 
been detected. This construction is implied whenever the 
original and the dual lattice are put into relation. The 
distribution of the values of this conditionally averaged 
Polyakov loop is shown in the different panels of Fig. |3| 
as histograms (with broad 12 bins, drawn in thick lines). 
For comparison, the distributions of the Polyakov loop 
over all lattice sites are shown (as fine-binned histograms, 
drawn in thin lines). One sees that, with /? approaching 
the phase transition from the confinement side, the dis- 
tribution over all lattice sites becomes flatter, whereas 
the distribution over the position of monopoles becomes 
increasingly shifted toward ±f . 

Within our set of runs above but not so close to the 
deconfinement temperature global flips of the Polyakov 
loop did not happen. Due to this effective Z{2) breaking 
the two plots for the deconfined phase show an unsym- 
metric distribution for the Polyakov loop over all lattice 
sites. At these temperatures the same asymmetry is seen 
in the distribution of the Polyakov loops for the monopole 
positions. 

In order to facilitate a qualitative comparison of the 
topological clusters with KvB caloron solutions with 
nontrivial holonomy L^o, we have defined an empirical 
asymptotic holonomy H for each lattice configuration. 
We have considered the sites on the lattice where the ab- 
solute value of topological charge density is less than the 
averaged absolute value of this density. We take the aver- 
age of the Polyakov loop over this set of asymptotic sites 
and call the resulting observable the asymptotic holon- 
omy H. The distribution of H for the five f3 values is 
shown in Fig. 0] 



^ For the purpose of this estimate the string tension was defined 
by the Polyakov loop correlator. 
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FIG. 3: (color online) The distribution of the Polyakov loop over sites where time-like Abelian monopole currents are detected 
(thick line). Monopoles are obtained by Abelian projection in MAG. For comparison, the distribution of Polyakov loops over 
all sites is shown (thin line). The upper row of figures shows the results for l3 — 2.2, 2.3, 2.4 (from left to right, confinement). 
The lower row of figures shows the results for /3 — 2.5, 2.6 (from left to right, deconfinement) . 



It can be seen from the different panels of this figure 
that the asymptotic holonomy H is concentrated near 
zero deep in the confined phase (/3 = 2.2 and 2.3) . 
The distribution is widened at (3 — 2.4 , already close 
to the deconfining transition, and it becomes again con- 
centrated around a maximum moving slowly between 
-1-0.5 and -f 1 for our temperature range in the deconfined 
phase. 

Note the sharp contrast characteristic for the confine- 
ment phase of the H distribution to the Polyakov loop 
distribution at the positions of time-like monopoles (as 
defined in MAG) in Fig. 13 In the deconfined phase 
the local Polyakov loop at the monopole positions is dis- 
tributed around the asymptotic (global) holonomy H. 

V. CLUSTER ANALYSIS OF SMEARED 
CONFIGURATIONS 

On one hand, for each of the four-dimensionally 
smeared configurations we have looked for clusters of 
topological charge. The topological density is assigned 
to the lattice sites according to the unimproved (naive) 
twisted plaquette definition. In order to form clusters for 
the subsequent analysis we have selected the sites where 
the absolute value of the topological charge density ex- 
ceeds some threshold value qc. The link-connected sites 



with q{x) > Qc or q(x) < —Qc form what we call posi- 
tive (negative) clusters of topological charge. Obviously, 
number and size of the clusters depend on the thresh- 
old qc- This value has been varied between the ensemble 
average of |g(a;)|, taken as lowest threshold, and a value 
taken 10 times larger taken as the highest threshold with 
the aim to get, for each smeared configuration, a maximal 
number of mutually disconnected clusters of topological 
charge. Therefore, this procedure finds an upper bound 
for the number of clusters per configuration (correspond- 
ing to the degree of smearing that we have chosen) . This 
number can be inferred from Table I and amounts to 20 
to 55 clusters per configuration, maximal at the lowest 
and highest temperature, minimal at T < Tc. Obviously, 
the number of clusters is not a well-defined quantity, and 
many of the shallow clusters should be ignored as mere 
extended background fluctuations. 

On the other hand, we have identified the complete 
monopole structure for each smeared configuration. This 
means that all monopole clusters have been found. A 
monopole cluster is a connected set of occupied links of 
the dual lattice, i.e. all links carrying non- vanishing mag- 
netic charge current. The density of monopole currents 
{i.e. the percentage of occupied links of the dual lat- 
tice) in the smeared configurations was about two orders 
of magnitude smaller than in the corresponding equilib- 
rium ones. Therefore, in the smeared ensemble monopole 
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FIG. 4: The distribution of the "asymptotic" holonomy H as defined in the text. The upper row of figures shows the results 
for (3 = 2.2, 2.3, 2.4 (from left to right, confinement). The lower row of figures shows the results for /3 = 2.5, 2.6 (from left to 
right, deconfinement). 



TABLE I: The numbers of topological clusters for 200 smeared 
configurations at different P identified as static dyons (D), 
non-static calorons (CAL) and without clear identification 
(OTHER). 



13 


T/Tdcc 


D 


CAL 


OTHER 


2.2 


0.54 


39 


624 


10323 


2.3 


0.68 


128 


755 


6512 


2.4 


0.94 


626 


511 


2983 


2.5 


1.32 


243 


84 


6989 


2.6 


1.87 


118 


20 


10204 



clusters were mainly non-selfintersecting loops, in the 
confinement phase yet percolating loops. From the set 
of all monopole clusters in a configuration, we have se- 
lected those isolated monopole loops that are either closed 
by periodicity in the time direction (then mainly formed 
by minimal length L = 6 static time-like monopole cur- 
rents) or closed within the four-dimensional volume (then 
mainly with the minimal loop sizes L = 4, 6, 8). If 
there was a complete covering of exactly one monopole 
loop (of either kind) with some topological cluster, we 
have classified this cluster either as a static dyon or as 
a nondissociated caloron, which are the two opposite ap- 
pearances of KvB calorons. Table I shows the corre- 



sponding numbers of topological clusters classified in one 
or the other way. Note that our classification, due to 
the very restrictive cut, has been passed only by a small 
fraction of topological clusters. The number of topolog- 
ical clusters in general is best defined closely below the 
phase transition (/3 = 2.4) where it reaches a minimum 
of w 15 per configurations. Approximately one third of 
them becomes successfully classified as dyons or calorons, 
the rest has no obvious interpretation in terms of the 
monopole content. The minimal monopole clusters men- 
tioned above could be intersecting with the large, perco- 
lating monopole loops, hence having escaped our atten- 
tion. 

With the 4-dimensional volume of 17.84 fm (at 
(3 = 2.4 the lattice spacing is obtained from aa^ = 
0.071(1) l^^) the density of identified calorons could 
be estimated as 0.143 fm"*, the density of dyons 
0.175 fm~^. The relative abundance extrapolated to the 
full number of clusters at this (3 would for the calorons 
result in 0.374 fm^'', still in the right ballpark given 
the constant topological susceptibility x ~ 1 fm^"* below 
Tdoc- Expressed in physical units, one can summarize the 
outcome as follows: the density of monopoles (dyons) 
continuously increases with the temperature. Near the 
phase transition the densities of calorons and monopole- 
like objects become comparable. We will discuss below 
to what extent the monopole- like objects arc really dyons 
and may be described as caloron constituents. The den- 
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sity of objects classified as calorons drops down above 
the phase transition. This does not come not unex- 
pected in view of the topological susceptibility decreas- 
ing with increasing temperature above Tdcc- It remains 
to be seen whether, in the deconfined phase, the config- 
urations with topological charge Q — ±1 really contain 
such a caloron. This would have to be expected from the 
scenario that calorons are mainly undissociated in the 
deconfined phase. 



TABLE II: The densities (in fm~*) of dyon and caloron clus- 
ters, as identified (in the second line extrapolated to the full 
number of topological charge clusters). 



u! according to Eq. as 



p 




D 


CAL 


2.2 


0.54 


1.14 • 10"^ 


0.0182 


2.2 


0.54 


0.018 


0.283 


2.3 


0.68 


9.78 ■ 10"^ 


0.0577 


2.3 


0.68 


0.062 


0.425 


2.4 


0.94 


0.175 


0.1429 


2.4 


0.94 


0.458 


0.374 


2.5 


1.32 


0.260 


0.0900 


2.5 


1.32 


5.55 


1.92 


2.6 


1.87 


0.514 


0.0872 


2.6 


1.87 


38 


6.45 



In order to corroborate (or relativize, for the decon- 
fined phase) this interpretation we have tried to charac- 
terize the monopole-tagged clusters by some cluster vari- 
ables. We have averaged the Polyakov loop inside the 
selected clusters again over all those sites where time-like 
Abelian monopole currents of cither sign are observed as 
part of the closed monopole loop that served the identi- 
fication. We call this quantity the Polyakov loop aver- 
aged over the "monopole skeleton" of the given cluster, 
(PL(Abelian monopoles)) cluster- From KvB calorons we 
expect this average for "static monopole clusters" to be 
close to ± 1 , whereas for "monopole-antimonopole pair 
clusters" it should be close to , the latter because of 
the internal dipole structure in terms of the Polyakov 
loop that reflects the positions of the constituents which 
otherwise would be hidden inside the single undissociated 
caloron cluster. We use this as the first cluster variable 
that should quantitatively characterize the type of clus- 
ter. 

As a second cluster variable we want to consider the 
topological charge assigned to the cluster. This total 
charge is difficult to assess because the topological profile 
extends below the threshold Qc- Applying this cutoff was, 
however, indispensable to localize the cluster. We have 
designed a model- dependent estimator for the total topo- 
logical charge of an identified topological cluster. This 
estimator works differently for the two types of clusters. 
For a static dyon we know from the analytic KvB caloron 
solution that its size depends on the holonomy parameter 
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(2) 



b is the inverse temperature, i.e. the period in time direc- 
tion. In the confinement phase assuming the holonomy 
parameter to be w = 0.25 the size of the dyons changes 
inversely to temperature. Actually, the sizes of the dyon 
clusters are distributed within some width. According 
to the caloron solution in the limit of separated dyons, 
the modulus of the topological charge density qmax in the 
center x^ax of the dyon cluster should be connected with 
the cluster size in the following way: 
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Ascribing by this relation a size to each cluster classified 
as static dyon we could obtain a cluster size distribution 
from the observed maxima of the topological charge den- 
sity inside the clusters. Then, considering such a cluster, 
we sum the actual topological charge over all sites x that 
belong to the cluster including the tail below Qc. That 
means that we sum the topological charge density q{x) 
within a tube with a 3- dimensional (spatial) radius R 
(related to f). This distance R should not be too large 
in order to avoid double counting of topological charge 
(by assigning sites to more than one cluster) and not 
too small (in order not to underestimate the topological 
charge in the tail of the cluster under consideration) . We 
use R — 3 f which, for the case of an isolated dyonic clus- 
ter (with an ideal, approximately exponential topological 
charge density profile), would estimate the total charge 
within 7% accuracy. This estimator corrects for the tail 
and serves to assign a topological charge to all clusters 
once they have been classified as static dyon clusters. 

An undissociated KvB caloron has a topological charge 
profile like that of an isolated ordinary instanton solution. 
In this case the maximum of the modulus of the topolog- 
ical charge density is related to the instanton size p as 
follows 



(4) 



Assuming that the clusters classified as undissociated 
calorons have such a charge profile we can obtain the in- 
stanton size p from the measured qmax of the cluster. 
Then we sum the actual topological charge density over 
all sites x inside a 4-dimensional ball with a radius 1.5p 
centered at x^ax- Finally, the result needs to be multi- 
plied by a correction factor 1.29 inferred from the exact 
instanton solution. In this way, we are in the position 
to define an estimated topological charge for any cluster 
once it has been identified as undissociated caloron. 

We show in Fig. scatter plots showing clusters of 
our topological clusters in the plane spanned by the es- 
timated topological charge of each cluster, Qduster ; 
and the Polyakov loop averaged over the monopole 
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FIG. 5: (color online) Scatter plots in the (Qduater, (PL(AbeUan monopoles))c;ustei ) plane for (3 = 2.2, 2.3, 2.4 (upper row 
from left to right, confinement) and for {3 = 2.5, 2.6 (lower row from left to right, deconfinement) . The circles represent dyon 
clusters, the triangles undissociated calorons. 
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FIG. 6: (color online) Scatter plots in the (Qduster, (P-L(Abelian monopoles))c!uster) plane for /3 = 2.6. The left figure is the 
sum of the other two. The figure in the center shows configurations (45 from 200 configurations) with topological charge = I 
(0.5 < IQI < 1.5), the right figure the complementary 155 configurations with zero topological charge Q = Q {\Q\ < 0.5). The 
meaning of the symbols is the same as in Fig. |S] 



skeleton, denoted as (PL(Abelian monopoles)) duster ■ that have been identified in one or the other way. 
The points in these scatter plots represent cluster As it can be seen from this figure, in the confined 
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phase with increasing /3 (here at /? = 2.3, 2.4) 
the two sorts of topological clusters are forming clus- 
ters on the scatter plot either closer to the points 
(PL(Abelian monopoles)) c!«ster = ±c, Qduster = ±1/2 
(dissociated) or (PL(Abelian monopoles)) duster — 
0, Qduster — ±1 (undissociated), with c ^ 1 with in- 
creasing temperature. In the confined phase at lower 
temperature, /3 = 2.2, clusters which contain closed 
monopole loops ( "undissociated calorons" ) are the most 
abundant objects (among the identified ones) and are 
forming a cluster, very similarly to /3 = 2.3, 2.4 . On 
the other hand, dissociated dyons are very rare and have 
not yet produced a pronounced pattern of non- vanishing 
Polyakov loop (when the Polyakov loop is averaged over 
the monopole skeleton of the cluster) . 

In the deconfined phase at /? = 2.5, 2.6 undissociated 
calorons are forming clusters in the same way as they do 
in the confined phase. However, they become increas- 
ingly rare objects with increasing temperature. Single 
(dissociated) dyons are presented now by asymmetric ob- 
jects: rare heavy dyons with an averaged Polyakov line of 
opposite sign compared to the holonomy H and frequent 
light dyons with averaged Polyakov line approximately 
equal to the value of the asymptotic holonomy. Undis- 
sociated calorons and heavy dyons are appearing almost 
exclusively in configurations with total topological charge 
equal to Q = ±1 while light dyons appear mainly in con- 
figurations with total topological charge equal to Q = 
and to a lesser extent in Q = ±1 configurations (see Fig. 
O. The charge values Q — 0, ±1 exhaust all cases among 
the smeared configurations on the given lattice size, that 
are found for the two /^-values in the deconfined phase. 
This suggests to view the topological background in the 
deconfined phase mainly as a gas of "spurious" (both 
with small action and topological charge) dyon-antidyon 
pairs. These pairs cannot be understood as forming to- 
gether a caloron. If they were related to calorons at all, 
more likely they have emerged from a caloron-anticaloron 
pair by annihilation of their heavy partners. For these 
dyons it has to be checked whether they are formed by 
locally selfdual or anti-selfdual field. Modeling the size of 
these objects starting from the topological density in the 
center is more questionable than in the other cases. Con- 
sidered only as static Abelian magnetic monopoles, they 
do not have to follow the dyonic interpretation. However 
they could be held responsible for the magnetic confine- 
ment. From this point of view it is no surprise that their 
density increases with temperature above Tdoc- 

Under further cooling, pairs of dyons and antidyons in 
this gas can eventually annihilate. In some cases (when 
the annihilation proceeds across the periodic boundary) 
this process ends up with Dirac sheets (i.e. constant 
magnetic fluxes) |2^ . These Dirac sheets are either stable 
or unstable depending on the value of holonomy |23| . 

With some reservations concerning the dyonic inter- 
pretation of the monopole-like clusters in the deconfined 
phase, we have made a comparison of the "dyon" radii 
found in the confined and deconfined phases. In the left 



half of Fig.[7|we show the size distributions derived from 
the gmaa; -values of clusters classified as single dyons. For 
the confinement phase, this yields a distribution with 
a single maximum as expected from the caloron model 
taken seriously. In the deconfined phase, the discrimina- 
tion between the abundant light dyon clusters with the 
same sign of the Polyakov loop at the monopole posi- 
tion and the rare heavy dyon clusters with opposite sign, 
always with respect to is reflected in the figure and 
highlighted by dotted and dashed lines, respectively. For 
comparison, we show in the right part of Fig.dthe p dis- 
tributions of isolated single caloron clusters found in the 
confinement phase. These distributions undoubtedly re- 
fiect the increasing size of the objects yet to be classified 
as single calorons. This tendency fits together with the 
increasing amount of dyons since, with increasing size, 
calorons eventually turn into discernible dyons. 

According to the relation (|2Jl, an estimate of w = 0.25 
for maximally nontrivial holonomy at /3 = 2.3, and 
uj — 0.1 corresponding to (H) — 0.8 at /3 = 2.6 could 
be taken for rough orientation. Since the inverse tem- 
perature is 6 = 6a (with — 6) in all our cases, we 
would estimate f w 2 x a(/3) for all dyonic constituents 
in the confined phase, whereas at /3 = 2.6 in the decon- 
fined phase the light dyons should have a typical radius 
f = 4.8 X a(2.6) and the heavy dyons (with lu replaced 
by tj = 0.5 — w = 0.4) a radius of f = 1.2xa(2.6). Notice 
that the lattice spacings for the two /3-values are widely 
different (they differ by a factor 2.75). Taking this into 
account the dyons at P — 2.3 (equally weighted in the 
confinement phase) would have a size of f « 0.6/ y/a, 
that happens to be almost equal to the size of the light 
dyons at /? = 2.6 (in the deconfined phase) which is 
f Ri 0.6A/y/a. The heavier dyons at f3 = 2.6 in the 
deconfined phase (which are actually rather rare) would 
be 4-times smaller in size, namely f ~ 0.16/y/a. For the 
lattice spacings a(/3) assumed here in order to express all 
sizes in terms of the SU{2) string tension at T = we 
refer to Ref. H. 

Concluding we can say that the distributions drawn 
in Fig.[7|do not agree quantitatively with the above esti- 
mates dictated by holonomy values H = L^o, but show at 
least qualitatively the expected pattern of size splitting. 
On one hand, the averaged holonomy H might not be a 
relevant parameter for the asymptotics of a single topo- 
logical cluster. On the other hand, we would not take 
these caloron-based estimates too serious for the static 
monopoles seen in the deconfined phase, because they 
are not observed in heavy-light pairs and are badly mod- 
elled by the formula for calorons and selfdual BPS dyons. 



VI. CONCLUSION 

We have analyzed the topological clusters of four- 
dimensionally smeared configurations at finite tempera- 
tures, across the deconfining phase transition, by means 
of studying the monopole-worldline cluster content of 
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FIG. 7: (color online) Left: distributions of the radii f of dyonic clusters for /3 = 2.3 (solid lines, confinement phase) and 
/3 = 2.6 (dashed/dotted lines, deconfinement pheise). The dashed line corresponds to heavy dyon clusters, whereas the dotted 
line refers to light dyons. Right: distributions of the scale-sizes p of isolated caloron clusters for /? = 2.2 (solid line), /3 = 2.3 
(dashed line), and /3 = 2.4 (dotted line). Both the values for r and p are presented in units of the lattice spacing a{/3). 



these clusters. This tool has allowed to identify part 
of the topological clusters either as static dyons or non- 
static calorons. We have to stress, however, that the 
number of separable topological clusters was maximized 
in the course of analysis such that a huge part of the clus- 
ters (more than 90 % at the lowest temperature and in 
the deconfined phase, approximately 70 % closely below 
the transition temperature) remained unidentified in this 
sense. This doesn't mean that they couldn't be classified 
using other means. The simple monopole clusters that 
were thought to identify clusters as dyons or calorons 
could have been intersecting with the percolating part 
that exists in the confinement phase. Another part of 
the topological clusters, mainly in the deconfined phase, 
should simply have been ignored because they are shal- 
low background fluctuations. These have nothing to do 
whatsoever with the well-localized, (anti)selfdual carriers 
of topological charge which had to be analyzed. 

On the basis of this separation we attempted a model- 
dependent reconstruction of the total topological charge 
for each individual cluster. The correlation between the 
Polyakov loop (averaged over the monopole skeleton) and 
the total topological cluster charge has revealed the fol- 
lowing pattern. At temperatures not too much below the 
deconfining phase transition, with the asymptotic holon- 
omy being still maximally nontrivial, H = 0, we see both 
nondissociated calorons and separate dyons, understood 
as the result of part of the calorons being already disso- 
ciated into dyons. The size of the nondissociated ones 
slightly increases with temperature. Dyons with differ- 
ent local values of the Polyakov loop are symmetric in 
mass and action in the confinement phase. The ratio of 
the number of dyons to the number of (yet undissoci- 
ated) calorons is increasing with temperature, such that 
we may speak about a process of dissociation with in- 
creasing temperature. At the transition temperature the 
densities of both objects has become comparable. 

Above the deconfining phase transition holonomy 
starts to become trivial, H ^ Q. Contrary to our ex- 



pectations topological excitations are not dominated by 
(undissociated or dissociated) calorons. Only half of the 
Q = ±1 configurations can be interpreted as undissoci- 
ated calorons. Judged solely from the average Polyakov 
loop at the static clusters, dyons are still prevailing above 
Tdec- The most abundant objects are light dyons and 
antidyons with static Abelian monopoles inside. That 
dyons are no more symmetric in mass and action would 
have been expected. But they do not appear anymore 
as pairs of light and heavy constituents in the same con- 
figuration. Therefore it is difficult to understand these 
configurations as dissociated calorons. 

The typical Q = Q configurations contains monopole- 
antimonopole pairs, each with only insignificant topolog- 
ical charge. Since the estimate of the charge rests on for- 
mulae for BPS dyons, some reservations about the topo- 
logical charge are due. Considered simply as static mag- 
netic monopoles (ignoring the topological charge cluster) 
together with other magnetic monopoles they could be 
held responsible for the spatial string tension (magnetic 
confinement) in the deconfined phase. From this point of 
view it is plausible that the density of monopoles does not 
follow the decline of the topological susceptibility above 

-^dcc ■ 

The logical next step will be to repeat the present 
study with a suitable fermionic topological charge den- 
sity applied to the (unsmeared) equilibrium gauge field 
configurations with the option of eliminating UV fiuctu- 
ations by mode truncation. 
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